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Abstract 

By recognising that stress-energy-momentum tensors are funda- 
mentally related to gravitation in spacetime it is argued that the classi- 
cal electromagnetic properties of a simple polarisable medium may be 
parameterised in terms of a constitutive tensor whose properties can in 
principle be determined by experiments in non-inertial (accelerating) 
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frames and in the presence of weak but variable gravitational fields. 
After establishing some geometric notation, discussion is given to ba- 
sic concepts of stress, energy and momentum in the vacuum where 
the useful notion of a drive form is introduced in order to associate 
the conservation of currents involving the flux of energy, momentum 
and angular momentum with spacetime isometrics. The definition 
of the stress-energy-momentum tensor is discussed with particular 
reference to its symmetry based on its role as a source of relativis- 
tic gravitation. General constitutive properties of material continua 
are formulated in terms of spacetime tensors including those that de- 
scribe magneto- electric phenomena in moving media. This leads to a 
formulation of a self-adjoint constitutive tensor describing, in general, 
inhomogeneous, anisotropic, magneto-electric bulk matter in arbitrary 
motion. The question of an invariant characterisation of intrinsically 
magneto-electric media is explored. An action principle is established 
to generate the phenomenological Maxwell system and the use of vari- 
ational derivatives to calculate stress-energy-momentum tensors is 
discussed in some detail. The relation of this result to tensors pro- 
posed by Abraham and others is discussed in the concluding section 
where the relevance of the whole approach to experiments on matter 
in non-inertial environments with variable gravitational and electro- 
magnetic fields is stressed. 

MSG codes: 83D05, 83G40, 83G35 

Keywords: Stress-energy-momentum tensor, constitutive relations, vari- 
ational, electromagnetic, polarisable, magneto-electric, gravitation, 
Maxwell's equations. 
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1 Introduction 



The laws of quantum-electrodynamics have been devised to describe the elec- 
tromagnetic interactions with matter according to the tenets of relativistic 
quantum field theory. However Maxwell's classical equations remain manda- 
tory for the description of a vast amount of natural phenomena. This ver- 
satility is in part due to the supplementary constitutive relations that are 
necessary to accommodate the wide range of materials that respond to elec- 
tromagnetic fields. Although in principle such relations can be derived from 
the underlying quantum description of matter, in many practical situations 
one must rely on experimental guidance to ascertain the classical response of 
materials to such fields. 

Once the unifying power of a spacetime formulation of physical phenom- 
ena became apparent with Einstein's relativistic world view, the natural 
mathematical tool for describing constitutive responses became the total 
stress-energy-momentum tensor for all matter and fields. Early suggestions 
by Minkowski and Abraham [2] for the structure of its electromagnetic 
component in simple media initiated a long debate involving both theoretical 
and experimental contributions that continues to the current time (see e.g. 
[3]),[1],[5],[6],[71,[8],[I0],[22],[23]). Although it is widely recognised that 
this controversy is an argument about definitions [TT] and that the relative 
merits of alternative definitions are undecidable without a complete (experi- 
mentally verifiable) covariant description of relativistic continuum mechanics 
for matter and fields, it remains important to clarify the many conflicting 
arguments that have appeared over the years and to offer new insights that 
may help in modelling the electromagnetic properties of moving media in the 
absence of a viable or complete description of field-particle interactions at a 
more fundamental level. 

Some way towards this goal is offered by (covariant) averaging meth- 
ods [18], [T7]. These however yield non-symmetric stress-energy- momentum 
tensors for electromagnetic fields in simple media. If the total stress-energy- 
momentum is to remain symmetric this implies that other asymmetric con- 
tributions must compensate and no guidance is offered to account for such 
material induced asymmetries. The need for a symmetric total stress-energy- 
momentum tensor is often attributed to conservation of total angular mo- 
mentum despite the fact that such global conservation laws may not exist in 
arbitrary gravitational fields. Although the magnitude of gravitational in- 
teractions may be totally insignificant compared with the scale of those due 
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to electromagnetism, gravity does have relevance in establishing the general 
framework (via the geometry of spacetime) for classical field theory and in 
particular this framework [13], [12] offers the most cogent means to define 
the total stress-energy-momentum tensor as the source of relativistic grav- 
itation. This in turn may be related to a variational formulation [16] of the 
fully coupled field system of equations that underpin the classical description 
of interacting matter in terms of tensor (and spinor) fields on spacetime. 

In this article stress-energy-momentum tensors are defined as variational 
derivatives and it is argued that the classical properties of a simple polaris- 
able medium may be parameterised in terms of a constitutive tensor whose 
properties can in principle be determined by experiments in non-inertial (ac- 
celerating) frames and in the presence of weak but variable gravitational 
fields. 

There has been a rapid development in recent years in the construction 
of "traps" for confining collective states of matter on scales intermediate be- 
tween macro- and micro- dimensions. Cold atoms and nano-structures offer 
many new avenues for technological development when coupled to probes by 
electromagnetic fields. The constitutive properties of such novel material will 
play an important role in this development. Space science is also progressing 
rapidly and can provide new laboratory environments with variable gravi- 
tation and controlled acceleration in which the properties of such states of 
matter may be explored. It will be shown below that the response of elec- 
tromagnetically polarisable media to such novel experimental environments 
offers a means to describe their electromagnetic constitutive properties and 
hence gain insight into the electromagnetic stresses induced by electro- 
magnetic fields in such media. Supplemented with additional data based 
on their mechanical and elasto-dynamic responses one thereby gains a more 
confident picture of the total phenomenological stress-energy-momentum for 
media than that based on previous ad-hoc choices. 

Throughout this article the formulation will be expressed in terms of 
tensor fields on spacetime with an arbitrary metric. Attention will be drawn 
to conservation laws when this metric admits particular symmetries. Thus 
the results have applicability to simple media in arbitrary gravitational fields 
and accommodate both media and observers with arbitrary velocities. 

After establishing some geometric notation, section 2 relates the electro- 
magnetic 1-forms e, b, d, h to the 2-forms F and G that enter into Maxwell's 
phenomenological covariant field equations in the presence of matter. Sec- 
tion 3 discusses stress, energy and momentum in the vacuum and introduces 
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the useful notion of a drive form that can be used to calculate electromag- 
netically induced currents involving the flux of electromagnetic energy, mo- 
mentum and angular momentum in Minkowski spacetime. In section 4 the 
definition of the stress-energy-momentum tensor is discussed with particular 
reference to its symmetry based on its role as a source of relativistic grav- 
itation. The constitutive properties of the media considered in this paper 
are delineated in section 5 in terms of a constitutive tensor on spacetime. 
This includes an account of general magneto- electric continua and leads in 
section 6 to a formulation of a self-adjoint constitutive tensor describing, 
in general, inhomogeneous, anisotropic, magneto-electric matter in arbitrary 
motion. The question of an invariant characterisation of magnto-electric me- 
dia is mentioned in section 7. In section 8 an action principle is established 
to generate the phenomenological Maxwell system and the use of variational 
derivatives to calculate stress-energy-momentum tensors is discussed in sec- 
tion 9. The computation of the electromagnetic stress-energy-momentum 
tensor, based on the action of section 8, is non-trivial for general media ex- 
hibiting anisotropy and magneto-electric properties in arbitrary motion and 
is presented in some detail. The relation of this result to tensors proposed 
by Abraham and others is discussed in the concluding section where the 
relevance of the whole approach to experiments on matter in non-inertial en- 
vironments with variable gravitational and electromagnetic fields is stressed. 

Notations follow standard conventions with spacetime modelled as a 4- 
dimensional, orientable, manifold M with a metric tensor field g of Lorentzian 
signature (— , +, -|-, -|-). TTM denotes the set of vector fields and FA^M 
the set of p— form fields on M. The set {e°, e^, e^, e^} denotes a local g- 
orthonormal coframe (a linearly independent collection of 1— forms) with 
dual frame {Xq, Xi, X2, X3}. If gab = g{Xa,Xb), the interior contraction 
operator ix^ with respect to Xa is written ia with = g°'Hx^, = gac^^ and 
summation over 0, 1, 2, 3. Metric duals with respect to g are written with a 
tilde so that X = g{X, -) G TA^M for X e TTM and 5 = g~^ia, -) G TTM 
for a G TA^M. The Hodge dual map associated with g is denoted ^^r. The 
following standard identities will be used repeatedly in subsequent sections 
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to simplify expressions. 
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2 Electromagnetic Fields 

Maxwell's equations for an electromagnetic field in an arbitrary medium can 
be written 

dF = and d G = j (8) 

where F G TJS^M is the Maxwell 2-form, G G TA^M is the excitation 2-form 
and j G FA^M is the 3-form electric current source 0. In general, the effects 
of gravitation and electromagnetism on matter are encoded in this system 
in -kG and j. This dependence may be non-linear and non-local. To close 
this system, " electromagnetic constitutive relations" relating G and j to F 
are necessary. In the following the medium will be considered as containing 
polarisable (both electrically and magnetically) matter with G restricted to 
a real point-wise linear function of F, thereby ignoring losses and spatial 
and temporal material dispersion in all frames. Continua endowed with such 
properties will be termed "simple" here. The electric 4-current j will be 
assumed to describe (free) electric charge and plays no role in subsequent 
discussions. 

The electric field e G FA^M and magnetic induction field b G FA^M as- 
sociated with F are defined with respect to an arbitrary unit future-pointing 
timelike 4— velocity vector field U G TTM by 

e = iijF and ch = iu-kF (9) 

"'^AU tensors in this article have dimensions constructed from the SI dimensions 
[M], [L], [T], [Q] where [Q] has the unit of the Coulomb in the MKS system. We adopt 
[g] ~ [L^], [G] ~ [j] = [Q], [F] = [Q]/ eo where the permittivity of free space eo has the 
dimensions [Q^ T'^M^^ L^^] and c denotes the speed of light in vacuo 
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Since g{U, U) = -1 

F = eAf/-*(cbAf/) (10) 

The field U may be used to describe an observer frame on spacetime and its 
integral curves model idealised observers. 

Likewise the displacement field d G TA^M and the magnetic field h e 
TA^M associated with G are defined with respect to U by 

d = iuG , and h/ c = % * G . (11) 

Thus 

G' = dAf7-*((h/c) Af7) (12) 

It will be assumed that a material medium has associated with it a future- 
pointing timelike unit vector field V which may be identified with the bulk 
4— velocity field of the medium in spacetime. Integral curves of V define the 
averaged world-lines of identifiable constituents of the medium. A comoving 
observer frame with A— velocity U will have U = V. 

3 Electromagnetic Stress, Energy and Mo- 
mentum in the Vacuum 

The historical development of Newtonian continuum mechanics led to the 
notion of a stress tensor in Euclidean 3— space that entered into the balance 
laws for momentum and angular momentum. With the advent of relativistic 
concepts this was generalised to a stress-energy-momentum tensor in space- 
time giving rise to conserved quantities in situations where the metric admits 
symmetries. 

The basic properties of the electromagnetic stress-energy-momentum ten- 
sor in the vacuumocan be succinctly discussed in terms of a set of "drive" 3— forms. 
In vacuo the Maxwell field system with a 3— form electric current source j 
satisfies 

dF = and eod ^ F = j. (13) 

^The notion of a classical vacuum here corresponds to spacetime devoid of all material 
{j = 0) although if j has compact support one can refer to "vacuum domains" where 
j = 0. All regions can admit non-zero electromagnetic and gravitational fields. 
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For any vector field Y on spacetime and any Maxwell solution F to this 
system one can introduce a "drive" 3-form associated with Y and F 

= ^(^^FA^F-zy* FAF) (14) 
2c 

This 3— form can be used to generate different types of conserved quanti- 
ties when the vector field Y generates (conformal) isometries on spacetimeU 
If /C is any (conformal) Killing vector on a domain of spacetime it then follows 
simply from the vacuum Maxwell-system above that 

dT^'' = -^z^FAj (15) 

Thus for each (conformal) Killing vector field these equations describe a 
"local conservation equation" ( drjc = 0) in a source-free region (j = 0). 
For K any unit timelike Killing vector one has from f|T^ 

r|^^ = A h A A' - 1{ eo g{e, e) + /io h)}^K (16) 

where h = /ig^b, b, e are defined with respect to U = K and /iq = 
The spatial 2— form e A h was identified by Poynting in a source-free region 
as proportional to the local field energy transmitted normally across unit 
area per second (field energy current or power) and |{ e^gie, e) + /iq g{h., h)} 
proportional to the local field energy density. More precisely J^^ tk is the field 
energy associated with the spacelike 3-chain S and J^a iRTK is the power fiux 
across an oriented spacelike 2-chain S"^. 

If X is a unit spacelike Killing vector generating spacelike translations 
along open integral curves then with the split: 

where iyQx = 0, the Maxwell stress 2-form nx may be used to identify 
mechanical forces produced by a fiow of field momentum current or pressure 
with momentum density 3-form Qx |19j . 

It is important to stress that different timelike Killing vectors give rise 
to physically distinct notions of conserved energy. For completeness the in- 
terpretation of "energy" requires further information related to its mode of 

■^i.e. In terms of the Lie derivative Cy i = for some scalar A. F is a Killing 
field when A = 0. Angular momentum currents follow in terms Killing vector fields that 
generate rotational diffcomorphisms. 
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detection. The existence of timelike parallel Killing vector fields (including 
those whose integral curves are geodesies) are further conditions that sin- 
gle out particular classes that may have priority in establishing appropriate 
notions of conserved energy. 

In general, in the absence of Killing vectors one loses strictly conserved 
currents (closed 3— forms) but a set of four local 3— form currents rf'*^ = 
can be defined in any local coframe. In any frame {Xa} with dual coframe 
{e''} the 16 functions T^*''^ = ixt * 'Ta^ may be used to construct the tensor 

usually referred as the stress- energy-momentum tensor associated with the 
above drive forms @ . 

The relationships between any stress-energy tensor T and the associated 
drive forms are given by 

Ta = MT{Xa,-)) and T = ^Ta A 6^)6" ® e' (18) 

In terms of the 3— forms Tc the symmetry condition T^c = Tcb is 

ec A Tfo = Cft A Tc (19) 



4 The total stress-energy-momentum tensor 

When spacetime contains domains with matter (where j may or may not 
be zero) such regions will in general have physical properties distinct from 
vacuum domains. 

If a coupled system of electromagnetic, gravitational and matter fields 
has a total stress-energy-momentum tensor 

rpTotal rpEM _j_ rpm+I (20) 

where T"^^^ describes matter and its interactions not included in T^'^ , then 
on general grounds, if T™'^' is symmetric, one has: 

V ■ T™ = (21) 

In view of the above comments on the role of particular timelike and spacelike Killing 
vectors in constructing conserved energy-power and momentum-force currents a more co- 
herent label for T might be the drive tensor 



7 



in terms of a (Koszcul) connection V on spacetime. Different authors parti- 
tion the total stress-energy tensor into a sum of partial stress-energy tensors 
in different ways. The divergences of certain partial stress-energy tensors are 
sometimes called pondermotive forces. 

If the connection V, induced from a connection on the bundle of linear 
frames over spacetime, is both metric compatible and torsion free, with grav- 
itational fields satisfying Einstein's equations, then T™"' must give rise to a 
symmetric stress-energy tensor T^°*'^' = T^°^°'K However any such symmet- 
ric tensor can be partitioned into non-symmetric partial tensors in infinitely 
many ways. Such a partition is then an expediency without fundamental 
significance. 

In theories of gravitation based on non-pseudo-Riemannian geometries 
the natural connection may have torsion. For example in an Einstein-Cartan 
theory with matter that gives rise to a connection with torsion, the gener- 
alised Einstein tensor Ein^'^ , determined by varying the generalised Einstein- 
Hilbert action with respect to orthonormal coframes, is non-symmetric and 
hence the source tensor T^'^' defined by 

Ein''^' = T^^ (22) 

is similarly non-symmetric. However for some forms of gravitational-matter 
couplings the variation of the total action with respect to the connection 
gives rise to algebraic equations for the connectionO In principle these can 
be solved for the connection which can always be decomposed into a sum 
containing the torsion-free metric-compatible (Levi-Civita) connection used 
in Einstein's pseudo-Riemannian description of gravitation. The generalised 
Einstein tensor Ein^^ can then be written Ein -|- 5" in terms of the Einstein 
tensor Ein and fl22l) becomes 

Ein = (23) 

where = T^^ — S is symmetric and divergenceless with respect to the 
Levi-Civita divergence. In such cases one may define the total stress-energy 
tensor as the source tensor for Einstein's equations fl23l) . It is then by defi- 
nition symmetric. If the natural connection V (determined by a connection 
variation of the total action) gives rise to dynamic torsion, determined by a 

^For example, locally SL{2, C) covariant couplings of spinor fields to gravitation fall 
into this category. 
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partial differential system involving all fields, the reduction to a geometrical 
formulation in terms of a metric and Levi-Civita connection becomes an im- 
practicality. In such a situation the definition of the stress-energy tensor is 
best left as T^'~''. This has two distinct divergences with respect to V since 
it is not symmetric. 

Such general considerations offer guidance in the construction of phe- 
nomenological partial stress-energy tensors based on either coarse-graining 
detailed interactions between fields or the introduction of effective degrees of 
freedom [13]. Indeed such phenomenological stress-energy tensors are often 
of greater value than actions based on "fundamental fields" since they can 
often be related more directly to experiment. Thus although in this article 
gravitation will be regarded as a background interaction the electromagnetic 
properties of a simple medium will be accommodated into certain constitu- 
tive tensors that respond to gravitation. We then demand that an action 
describing such a medium in the absence of free charges give rise by vari- 
ation to Maxwell's phenomenological equations for a simple medium and a 
symmetric stress-energy tensor. 



5 The constitutive tensor for simple media 

In general G may be a functional of F and properties of the medium^. 

G = Z[F,...] (24) 

Such a functional induces, in general, non-linear and non-local relations be- 
tween d, h and e, b. These relations may be explored either empirically or 
by coarse graining a suitable macroscopic model. For general linear continua 
one may have for some positive integer and collection of constitutive tensor 
fields Z on spacetime the relation 

G = Sf^QZW(V"F,...) (25) 

in terms of some spacetime connection V. Additional arguments refer to 
variables independent of F and its derivatives. In this article, for the simple 
linear media under consideration, we restrict to 

G = Z{F) 

^ e.g. electrostriction and magnetostriction arise from the dependence of Z on the 
elastic deformation tensor of the medium |15j . 
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for some constitutive tensor field Z . In the vacuum G = e^F. 

A particularly simple linear isotropic medium may be described by a bulk 
4— velocity field V , a relative permittivity scalar field e and a non-vanishing 
relative permeability scalar field /i. In this case Z follows from 

G ~ ~ 

— =e iyF A 1/ - //-^ * {iv^F A V") (26) 
Co 

=ie--)ivFAV--F (27) 
^ /i 

In a comoving frame with U = V ^I7\ becomes 

d = eoee and h = (yUo/i)~^b (28) 

For a non-magneto-electric but anisotropic medium, the relative permit- 
tivity e and inverse relative permeability become spatial tensor fields on 
spacetime. Thus e : VK^M VK^M and /x"^ : TA^M ^ TA^M for all 
a e TA^M where 

e(i7) = 0, 2ye(a) = 0, /i-i(i/) = and iv^l~^{a) = Q. (29) 
The more general constitutive relation is then given by 

— = e{ivF) ^V - ^(/i-^(zy*F) A V) (30) 

^0 

which in the comoving frame with U = V becomes 

d= eo e(e) and /ioh = /i~"^(b) . (31) 

Based on standard thermodynamic arguments the inverse relative permeabil- 
ity and relative permittivity tensors are symmetric with respect to the metric 
9- 

iae{P) = ij^e{a) and isyU~"^(/5) = z^/i~^(a) for a,PETA^M (32) 

In general, the electromagnetic fields may be related by 

d =C'%e) + r (b) 
h =C'^(e) + C'^(b) 
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where C'^^ C^^, C^^ : TA^M TA^M are spatial tensors satisfying: 

C(V^)=^v^(C(«))=0 (34) 

and therefore 

Cinvia)) = nviCia)) = ({a) (35) 

for ( = C^'^, C^^, and for all a G TA^M, where vry projects spacetime 

1-forms to spatial 1-forms with respect to V, on spacetime 

TTv : TA^M TA^M , wy = ldi + V ®V (36) 

From fl3Tl) . fl33l) it follows that if = C^*^ = in some frame then ('^'^ = Cq ^ 
and (^^ = (/io/i)~^ in that frame. For such materials,however, one cannot 
assert that ^'^^ remain zero in all frames. Media with constitutive relation 
( !33|) are often referred to as magneto- electric [20]. We prefer to use this term 
to describe intrinsic magneto-electric media and will return to this point in 
section [71 

The tensor fields C^^, and (^^ are encoded into the tensor Z : 
TA^M TA^M such that G = Z{F). Since 

Z{a + (3) = Z{a) + Z{(3) and Z{Xa) = XZ{a) (37) 

for all A e TA^M and G TA^M the constitutive relation may be ex- 
panded in a local co-frame field {e°, e^, e^, e^} as 

iGabc" Ae' = IZ'^bF^ae" A (38) 

where 

vcrf lycd lydc rydc 

^ ab — ha — ah — ^ ha 

These conditions alone imply that the tensor Z has 36 independent com- 
ponents, although additional symmetry conditions given below will reduce 
these to 21. From the definition of G in terms of comoving fields and (l33l) . 
the relationship between Z and {C^'^j C'^''; C'''^; C'^''} follows as 



Z{F) =('^%iyF) AV + C^^{iv*F) A V 

- HC^'iivF) AV)- ^(C"^(zv ^F)AV) 
and hence by contraction with V 

C'^O = WZ{^ A Vh C''(0 = -^vZ{^i^ A V)) ,^ 
C^<'(0 = ^ A V) , C''(0 = -V ^ Z{^{^ A V)) 



(40) 



(41) 
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6 Symmetry of the constitutive tensor. 

The adjoint of any tensor T : TA^M TA^M, is the tensor : TA^M 
TAPM defined by 

aAi.T{l3)=(3Ai.T\a) ioTa,(3eTAPM (42) 

Clearly T^t = T. If p = 1, (gZ]) gives 

isT{(3) = i^T\a) for a,(3e TA^M (43) 

The symmetry conditions for the relative permittivity and inverse perme- 
ability tensors imply that ('^'^ and (^^ are self adjoint. This symmetry is 
generalised to magneto-electric media: 

^dct _ ^dc ^ ^hbt _ ^hb ^dbt _ _^he ^44^ 

i.e. Z is assumed self-adjoint 

Z = Z^ (45) 

or, raising indices with the metric 

^abcd ^ ^cdaft ^^Q^ 

Using sequentially (|4T]),(ll,(ll5]),([l,([6]),(|iI]),([I]),([l this condition 
yields 

a A ^C^^{(3) = - a A MyZ^p AV)) = aAV A *Z(*(/5 A V)) 

=-k{(3 AV) A -kZ{a AV) = iy^p A ^Z{a AV) = i^ii A iv*Z{a A V) 
= ★ /3 A C^'(a) = A = -/5 A ^C^%a) = -a A <^"V) 

i.e. = — C^'^^- The remaining equations in (l44l) follow similarly. 

It follows from (!39|) and (l46l) that the number of independent components 
of Z reduce from 36 to 21. 

7 Intrinsic magneto-electric media and the 
Post constraint. 

A constitutive tensor Z describes a non intrinsic-magneto- electric medium 
if there exists a velocity field V for the medium such that = and 
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( = 0. Thus a constitutive tensor Z is intrinsically magneto-electric if 
there does not exists a velocity field V such that = and (^'^ = 0. If 
Z{F) is decomposed with respect to an arbitrary frame U ^ V one may 
find all tensors ( non-zero, even for media that are not intrinsically magneto- 
electric. For a general constitutive tensor, it is a matter of linear algebra to 
decide whether it describes an intrinsically magneto-electric medium or not. 

A useful characterisation of magneto-electric media may be given in terms 
of invariants constructed from Z and the metric. One such invariant intro- 
duced by Post [21], [9], is 

X = ta^b*iZ{e^Ae')) (47) 
In terms of spatial tensors with respect the the medium velocity V 

X=^a^b*{Z{e'' Ae')) 
=^a^b^{C''''{^v{e'' A e')) AV)+ laib * (C*(^y*(e" A e')) AV) 

+ ^a^b{C'''Me'' A e^)) AV) + ^a^b(C'^^v^*(e'^ A e")) A V) ^^^^ 
=Mb * (C''(^y^(e" A e'')) AV) + ^a^b{C'\^v{e'' A e'')) A V) 

since 

^a^AC'"'Me'' A e")) A V^) = latbiC'iiv^e^ A e^)) A V') = 

Using 

^(e" Ae^ AV)= ^(e" A A e")K = K * (e" A A e^^ A e'^)ed = -e'^^'^V^ed 
and 

A ■(/ A efe A ta) = A e;, A ^ A = *(ea A Cf, A A ef)V''if = -SabefV^i 
with Sabef^"'^'^'^ = S'^Sj — 6^5^, the first term on the last line of fjiHl) yields 

iJs * iC'^^iivMe'' A e')) AV) =i. (^^''(^(e" A A i/)) A A e, A ej 

= -i.{aer)AV AesAea)e'''''^V, 
=i^C''\er)easefe'''''V''V = 2zaC''(e'^) 

while the second term is 

M6(C"'(^y(e" A e^)) AV) = zMC^'^v'^e') AV) - v,(C'^°(^;V) A V) 

= 2zv^z,(C^(e'') AV) = -2z,C''(e^) 
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Hence using (jUj) 

X = 42.C''(e'^) = -42aC'''(e") (49) 

Thus since V is an arbitrary medium velocity, a sufficient condition for a 
medium to be intrinsically magneto-electric is that x 7^ 0- 

However some intrinsically magneto-electric media may have x = 0- Fo'^ 
example, consider the self-adjoint constitutive tensor given, in some local 
orthonormal coframe {e°, e^, e^, e"^}, by 

Z{F) = ^236° A + Fi3e° A - Foae^ A - FqiC^ A 

Then with V = Xq 

r (0 = (^lOe' - {^20e' and C'^(0 = -{H^e' + (^sOe' 

and X = 0. However one easily verifies that ^ with respect to any 
arbitrary unit timelike V . Hence Z describes an intrinsically magneto-electric 
medium. 

A minimal set of invariants whose non-vanishing is a necessary and suffi- 
cient condition for a medium to be intrinsically magneto-electric is not known 
to the authors. 

8 Action for Source Free Electromagnetic Fields 
in a Simple Medium 

The classical equations describing the total system of matter and fields will 
be considered as arising from the extremum of some total action functional 
under suitable variations with compact support. This action should be con- 
structed from an action density 4-form on spacetime in terms of (pull-backs) 
of sections (and their derivatives) of field bundles carrying representations 
of local symmetry (gauge) groups and maps between them. Observed local 
symmetries in nature arise in such a formalism by ensuring that the action 
4— form is a scalar under local changes of section. To maintain these covari- 
ances appropriate connections are required to define tensorial (and spino- 
rial) covariant derivatives of sections. In addition the action may depend on 
tensor-valued functions of these sections. All variational principles require a 
specification of what objects in the action are to be varied and these then 
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constitute the dynamical variables of the theory. In the following we con- 
centrate on a contribution A to the total action arising from the effects of 
the electromagnetic field and gravitation in different types of "media". We 
exclude from this A the interaction with charged matter and the dynamics 
of the gravitational field itself. Included is the effect of the electromagnetic 
field on a polarisable and magnetisable medium assumed to be described in 
terms of a particular constitutive tensor Z . In particular we explore how 
the response of the medium to gravitation as well as the electromagnetic 
field can be used to establish the stress-energy-momentum tensor associated 
with different choices of constitutive tensor. Thus the action 4-form A will 
be taken to depend only on the spacetime metric and the class of Maxwell 
1-form potentials A with F = dA. The dependence of the tensor field Z on 
these variables will be explored in some detail below. 

We have insisted that in the absence of free charge the electromagnetic 
fields F and G for a simple medium in any spacetime metric must satisfy 

rfF = and rf*G = (50) 

Before generating an electromagnetic stress-energy- momentum tensor from a 
particular contribution to the total action it is necessary to verify that these 
field equations arise by suitable variation. Consider then the contribution 
S[A, g] = Jj^j A where F = dA,G = Z{F) with Z = Z^ and 

cA = iF A^G = iF A*Z(F) (51) 

If a prime denotes the variation with respect to A, then working modulo d: 

cA' =1 i^dA' A i.Z{dA) + dAA i.Z{dA')^ 

=dA' A i<Z{dA) = A' Adi< Z{dA) =A'Adi<G 

Hence the source-free Maxwell equations (|50l) follow by variation with respect 
to A from the action ( 15T1) . Note that the symmetry condition (H5il of the 
tensor Z is essential in this variation. 

9 Variational derivatives and stress-energy- 
momentum tensors 

To effect the metric variations of the above action functional let t ^ gt he 
a curve in the space of Lorentzian signatured metrics, with go = gt\t=o- The 
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"tangent" to the the curve t ^ gt aX the point t = is written g: 

gt = go + tg + 0{f) (52) 

For a general object K which may be a tensor or a map which depends on 
the metric g, write similarly t — > Kt as the one parameter set of objects 
encoding the dependence of K on gt, Kq = Kt\t=o and K = -^Kt\t=o, so 

Kt = K^ + tk + 0{f) (53) 

Kt will be referred to as the metric induced lift of K. 

One may represent the local variation gt in different ways. One way is to 
vary the components of gt with respect to a fixed local co-frame {cq}, i.e. 

gt = {9t)abeo ® Co where {gt)ab = ^t((^o)a, (^o)b) (54) 

One can set the fixed frame to be orthonormal with respect to the unvaried 
metric so that {go)ab = Vat = diag(— 1, +1, +1, +1). The derivative g is 
therefore given by 

9 = Qabel ® (55) 

Alternatively one may vary the co-frame: i.e. choose a one parameter set of 
coframes t — > for a = 0, .., 3 such that |t=o = Cq and 

gt = Vabct ® e'l (56) 

The derivative of t ^ at t = follows from 

e^ = e[; + te" + 0(t2) (57) 

The derivative g may also therefore be written 

9 = Vab (e" (S)el + e^Q® e') (58) 

The drive 3-forms Ta associated with any action 4-form A are defined by the 
variation of A with respect to the orthonormal coframe as 

A = A r„ (59) 
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If the variation of A with respect to the ortho-normal coframe is induced 
entirely from the metric g (and the metric compatible torsion-free Levi-Civita 
connection) then 

Ta = 2lx, ( 1^ ) r/ac (60) 



This follows immediately by equating ( l55i) and ( l58ll : 

Qabel ®el = r]ab (e" ® + ® e^) 

so 

9ab = Vcd {e%Xa)6^ + S^a^'^iXh)) = ea{Xk) + eb{Xa) = IbCa + ^aCf, 

since one may drop the subscript here without ambiguity: Xa = {Xo)a- 
Then 

A =e° ATa = ^9ab = ^(ea(Xfe) + eb{Xa)) = 2-^ea{Xb) 

OQab Ogab Ogab 



--2ea A ib 



6A 



.^9ab, 

By (ITSl) the tensor associated to the is given by 



-2^[-^\ea®eb (61) 



Jab 



and is manifestly symmetric. 

In the following it is necessary to make explicit the metric dependence 
of various elements that enter in the action 4— form A and in particular to 
pass between vector fields and forms using the varied metric gt- Thus the 
notations gt : VTM TA^M, X ^ gt{X) and g^^ : VK^M TTM, 
a I— > g^^{a) for the metric dual of vectors and 1-forms with respect to the 
metric gt are used. For vectors or 1-forms which already have a subscript 
or t we continue to use the tilde notation without ambiguity so that, for 
example, the gt metric dual of the vector Xf can be written Xf = gt{Xt). 

Following flS51) one has the maps *t a-iid ^ and from the Leibnitz rule 
(evaluated at t = 0) 

(i,a)' = lira + *d (62) 
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for all at G Y!S?M . It follows simply (see appendix) that 

= A * a - ^(6" A iaO) for at G TA^M (63) 
Taking the derivative of $ A = \1/ A with respect to t gives 

$ A = ^ A for $, ^ G TA^M (64) 

Thus with the metric induced lift of the constitutive tensor Z: 

Zt = Zo + tZ + 0{t^) (65) 

one writes: 

{^Zt{F))- = ^Z{F) + ^Z{F) (66) 

Since there is a one parameter set of Hodge duals, we need to distinguish 
and to- Furthermore P2|) becomes 

a A *tT(/3) = /5 A *tT^' (a) (67) 

for all a,f3 e TA^ M and (US]) becomes 

for all a,/? G rA^M. 

10 Computation of the stress-energy-momentum 
tensor 

In this section the variation of the above action (ISTl )is explored for a partic- 
ular choice of the metric dependence for Zt, corresponding to a perturbative 
response of the medium to gravitation. 

For a general lift the action 4— form (IHT]) is written 

cAt = lF AMZt{F) (69) 

hence 

cA = i(F A^Z(F) + F A^Z(F)) (70) 
18 



From 0631) 

F A =F Ai.G = F Ae" Aiai^G - F A *(e" A laG) 

=e^ A F A iai^ G - A iaG A *F (71) 
=6" A {F Aia^G - iaG A -kF) 

and 

F A i<Z{F) = 2FA i.^{F)ibe'' = 2e'^ A z J F A ^t— (F) 
Therefore the drive forms are given by: 

CTa = lF Aia^G- \iaG A^F^iAf A *^(F) ) (72) 

For a physical medium with bulk motion that can sustain elastic stresses 
associated with its atomic constituents one expects that the history of such 
bulk motion should have some influence on the constitutive properties via 
some associated 4-velocity fielcS. To include the possible dependence of the 
stress-energy-momentum tensor on such bulk motion of the medium one re- 
quires Z to depend on this motion in some manner. In fl33p Z is specified 
in terms of electromagnetic fields measured in the comoving frame V of the 
medium. It is therefore natural to prescribe a lift of this expression involving 
the lifts of Vo and {Co'^, Co*^, Co°5 Co^}- The natural lift of the medium velocity 
Vq is 

Vt = , (73) 

The metric dual of Vt is given by \4 = 9t^t) and the projection Try ( |36!) is 
lifted to 

VTi = Id4 + yt®Ft (74) 

The decomposition ( 14T|) of Zq and Zt with respect to the medium veloc- 
ities H and Vt is given by {Co^Co^Co^C} and {Cf , Cf , C^^^ } respec- 
tively, following the notation (!53|) . 

The lifted tensors C^^*', C,f° ^ Q'^ and will be now chosen to satisfy three 
properties: 



^Relativistic strings and membranes with dynamics that arise from re-parameterisation 
independent actions are an exception since, without "constituents" , no preferred 
parametrisation of their histories should be identified. 
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For all t in the neighbourhood of t = 

C,|*=o = Co for a = Cf,Cf,Cr,Cf (75) 

For all t in the neighbourhood of t = they map the vector space that 
is (74— orthogonal to Vt to itself. This is achieved by lifting ([3i 



aVt) = and tvA{a) = (76) 

for Ct = ,Cf^,C^,C^ and all a e TA^M. The corresponding lift of 
§B is 

C(vr,(a))=vr,(C,(«)) = C(a) (77) 
For all t in the neighbourhood of t = they retain the adjoint conditions 



(Cf)^'=Cf% (Cf)^'=C'' and (Cf)t* = -C^ (78) 
These requirements are all satisfied by setting 

^xCt(«) = I {^xCoi'^ta) + Zg-.^{C,y"{7TtgtX)'^ (79) 

for all a G TA^M and X e FTM, i.ejfl 

^xCf^(«) =1 (^xC(vr*«) +^^-i„C(vr,<7,X) 
^xC''(«) =1 (^xC(vri«) +^,-i,C(vr*^7*X) 

ZxCf («) =1 (^xC(^*«) -V^-^o^(^*^*^)) 



(80) 



* For an isotropic, non-magneto-electric medium ([27]) and ([28]) . the lifts (|80|) reduces 
to the lifts 

Cf = 67r*, C^ = ^-\t, Cf = and = 
which in a comoving frame yield the relations 

dt=eoeet and ht (^o M)~^bt 
where the scalars e and /i^^ are independent of the ambient metric. 
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To verify that 0791) obeys fl75p note that from fl36p 
go'nogoX = X + go{X, Vo)Vo 
and hence from fl68l) and fl34l) 

Thus at t = ([79D becomes 

ixCt{(y)\t=o =\ (^zxCo(7ro«) +^g-i«(Co)^"(7ro^o^)) = «xCo(a) 
using ( l35|) . 

To verify that ( 1791) obeys ( 1761) observe that 

Now 7rt\4 = so the first term vanishes. Also using (l68l) 

=a/-5*(K„K,) V^-...xCo(K,) = 
Hence Ct(^) = 0- Likewise 

^y.C(«) = I (v.Co(vr<«) + V^JCo)^°(vr*^<)) = 
Finally to verify that ([80]) obeys ([78]) use ([68]) and ([80]) twice 

In a similar way it follows that (Cf = Cf' and (Cf^^)^* = Cf" . Thus ([HO]) 
provide natural conditions for the lifts ( 1751) to (178li^l. 

^ The requirements (|75ll78p are not meant to be exhaustive. Other lifts could involve 
gradients of the spacetime metric corresponding to gravitational tidal effects on the con- 
stitutive tensor. For example if 7?. is the curvature scalar associated with the Levi-Civita 
connection then the lifts 

also satisfy (1751178^ . 
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Inserting the relations fl80l) into Z{F), fl4UI) . the action 4— form (1^^ be- 
comes 

2 c A, =F A MCf'itv.F) AVt)+FA ^,(Cf (2y,^*F) A Vt) 

To ease the density of notation in the following, the symbol now stands 
for cbf and stands for ^. The lifts 

et = ivtF, ht = ivtMF , et = gt^^iy^F) and = g~^{iv,-ktF) (82) 
satisfy 

TTtet = et and vr^b^ = ht (83) 

Sequentially using (13]),(IS]),(Il,(IHn]),(IM]),(ll,(ISl),(I3]),the first term on the 
right hand side of 0811) becomes 

F A MCti^t) A Vt) =F A tvMtiet) = -et A ^tCti^t) = -(*a)^i.Cf^(e*) 

= - K^il) {t~.Xti^tet) + ^~.Xti^tet)) = -{^tl)^~.Xti^t) 
= -etA MCti^t) = FA tvMo'i^t) = FA MCti^t) A Vt) 

Similarly sequentially using (Il,(l6]),(|4]),(IHQl),(IHl,(!l,(l6]),(l3^ 
the second term on the right hand side of ( IHTl) gives 

FA^tiCt\^t)AVt) 

=F A tvMt\ht) = -et A ^tCti^t) = -i^tl)t~.Xt\ht) 

= - IKI) {^~.Xo\^tht) - ^^So''i^tet)) = -^mI) (^s.C(b*) - ^5,Co'^(eO) 

= - |e, A MCo\ht) + Iht A ^,Co°(eO = ^F A ^vMo\^t) - \mF A ivM^^^t) 

= \F A Md^t) A Vt) - \M{C:^\et) A V) A mF 

= \F A *i(C(bt) A Vt) -\FA MMiCi^t) A V) 

= \F A MiCi^t) A Vt) + \FA C^'iet) A V 

It is useful to record from this calculation that 

-U^tl) {^~.Xtiht) - ^iXo''i^t)) = \FA MiCi^t) A V) + \F A Ci^t) A V 

(84) 
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Sequentially using ©, O, ©, (HD, dHH]), dM]), dBl the third term on the 
right hand side of (ISTI) yields 

FAC^(eJ 

= - C^(ei) ^Vt^ MMF = C^(eO A Mht = A ^iC^(eO 

=i(^^l)(^b,Co'^(e^)-^i.C(b^)) 

= iF A ^i(Co*(bi) A Vt) + iF A C^(e,) A 

Finally on sequential use of p. (|2l). (HI). (|80|). (|83|) 

F A cf (bt) ^Vt = - Ctiht) ^Vt^ mmF = (b*) a ^^b^ = b^ a ^^Cf (b*) 

= Kl)%C,''(b,) = (^b.C(^*b,) +zg^C(vr,b,)) 

= Kl)^b.Co''(bO 

and so by reversing this sequence of steps 

F A (bi) A \4 = F A C(bi) A \4 
Hence ( ISTi) simplifies to 

2 c =F A A y,) + F A *,(Co*(^y,^tF) A V,) 

+ F A Co^''(zv^F) A \4 + F A C(v,*ti^) A \4 

i.e. the constitutive tensors C,t in the action may be replaced by Co and hence 
the metric dependence of is seen to reside solely in Vt and Vt- 
The derivative of (!85l) at t = is given by 

c A =F A ^(C'^'=(zvF) A + F A i.{C^\iyF) AV) + F A *(C'^°(zyF) A V) 

+ F A ^(C'^'^(2y*F) A + F A ^(('^^(v * F) A y) + F A *(C'^'^(iy^F) A 
+ F A ^(^^^(iy^F) A + F A C'(vF) A ^ + F A C'(iyF) A V 
+ F A C''^(V * F) A r + F A C''^(^y^i^) A V + F A C^^(iy * F) A 

(86) 

where the subscript is omitted on the right hand side. 

To determine the drive forms, observe that there are three different types 

of term in flS^ which contain i<, V or V . Since V = {g{V))' = g{V) + g{V), 
terms in ic, V and g{V) can be collected to give: 

2cA = Ai + A^ + A3(y) (87) 
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where 

2 c Ai =F A ^C^^ivF) AV) + FA ^C^^{ivi^F) A V) 

+ F A i.{C^^{ivi^F) AV)+F A ("^^{iv^F) A V , 

2cAy=FA i.iC'^'iiyF) AV) + FA *(C^°(2yF) A g{V)) 

+ FA ^(C*(V ^F)AV) + FA MC^\tv^F) A giV)) 

+ F A C (V^) A y + F A C (^yi^) A giy) 

+ FA C^^(v * F) A y + F A C''^(iy * F) A g{V) 

and 

2 c Ag^v) =F A ^{(""{tyF) A g{V)) + F A ^(C*(^y^F) A ^(V^)) 
+ F A C^^(iyF) A ^(V) + F A C^^(iy * F) A ^(r) 

The third term on the right hand side of fl88p may be expressed as 

F A ^(C'^^(iy^F) A V) =F A iy ★ (C'^'^(2y^F)) = -iyF A <*(iy^F) 

=iy^F A <^"(2yF) = -^F A iy ^ C^"(^y^) 

= - iy ★ C'''(vi^) A iF = -F A ^ * (C'''(^\/i^) A V) 

using sequentially ©, ([6]), dH]), ([6]), ([I]), ([61]). The fourth term on the right 
hand side of flHHl) may be expressed as 

F A C^^iiv^F) AV = -F AV A (^^{iyi^F) = -F A V A ^ ^ C'^^{^v^F) 
= - i^C^^iiv^F) A ^(F AV) = iv^F A ^C^^^y^i^) 
=2y^F A <^^(iy*F) = -^F A iy ^ C'^iiv^F) 
= -iv^ C^'^iiv^F) A ^F = -F A i * (C^'^iy^F) A 

using sequentially (P, ©, ([Ml), ©, dH, ([6]), ([1]), ([611). 
Hence from (17T1) 

2 c Ai =F A ^(c'^'=(iyF) A i/ + ('^^(iy^F) A V - ^((''^(^vi^) A i/) - *(C^^(iy*F) A V)) 
=F AiG 

=e'' A (F A ia * G - A -kF) 

(91) 
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To collect terms in Ky observe that by differentiating (175]) . V = XV where 
A = e^{V) Va one has 

2 c Ay =2A(F A^(C^"(zyF) A V) + F A ^C^^izy^F) A V) 

+ F A C^"(iyF) AV + FA C^^(iy ^F)AV) 
=2XF A = 2zye" Ki^ A = 26"^ A A 

The first two terms on the right hand side of fl90p become 

F A MC'VvF) A giV)) + FA ^iC'^^v^F) A giV)) 

=F A ^{ivG A g{V)) = ivG A g{V) A*F = -g{V) A lyG A ^F 

and the last two terms on the right hand side of (!90|) become 

F A C^'iivF) A g{V) + F A C^^(iy ^ F) A g{V) = -g{V) AF Aiy^G 

so using g{V) = e"-{V) + e" ea{V) = 2e" K + ivi^a A e") one has 

2cAg(v) = - g{V) A {ivG A^F + F Aiv^G) 

= - (2e" Va + iv{ea A e'*)) A {iyG Ai<F + FAiv^G) 

= -2e"' A {lyG Ai^F + F Aiv *G) + Ca Ae"" A iv{ivG A-kF + F Aiy -k G) 
=e A ( - 2{ivG A-kF + F Aiv -kG) + e'' A {iyF Aiy^G- iyG A ^F)) 

Adding this to 2Ay gives 

2cAy + 2cAgiy) =e A {2Va{ivF Ai^G- iyG A i^F) + e'' A e" A {iyF Aiy^G- iyG A i^F)) 
=2V;e" A (iyF A - iyG A ^F) - e'^ A A iviiyF A *G - iyG A 

Using the relation -kG = iv'kGAV + iv-*^ iyG and the similar relation for 
•kF, the combination above may be writen 

iyF A^G- ivG A -kF 

=ivF A iv^G AV + iyF A iy^iyG - iyG A iy^F AV - iyG A iy^iyF 
= ★ s — iviiyF A ^iyG — iyG A MyF) = -kS 

where the 1-form 

s = ★ fiyF A iv-kG AV + iv-kF A iyG A V^ (92) 
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Hence 

2cky + 2c kgiv) =6" A (2K*s - A iy^s) (93) 
Adding together (19T1) and (l93l) gives finally: 

2ck = e ^{F Mai^G -iaG ^i.F + 2K*s - e'* A (94) 
Hence the drive forms are given by 

CTa =\{F Mai^ G - laG A i^F) + K 7»r S - A * » (95) 

with associated stress-energy-momentum tensor: 

T=\ [i^F ® i^G + iaG®z^F -i.{F A i.G)g + V ^s + s® (96) 

The tensor T above coincides in Minkowski spacetime with that attributed 
historically to Abraham. It is derived here in a considerably wider context. 
In terms of comoving fields the drive forms can be written: 



(97) 



(98) 



c Ta =Vaiv (e A T^d + h A ^^b) — ^ (e A ia ^ d + i^d -k e) 

- i(h A * b + i^b * h) + 2vae AhAV + CaAeAh 

and hence 

T = - |(e O d + d ® e) - |(h ® b + b O h) 

+ l{g{e, d) + g{h, h)){g + 2V ^V) + {V ^ S + S ®V) 

where the Poynting 1-form 

S = i.{V AeAh) 

One may express the expressions above in terms of comoving polarisation 
1— forms p and magnetisation 1— forms m, defined in terms of comoving 
electromagnetic fields by 

d = e + p (99) 
h = b-m (100) 

Thus 

G = F + V (101) 
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where 

p = pAl/ + *(mAl/) (102) 

Then one finds 

T-c = rl + tI + tI + 

where 

Icrl =i^i.G NF - i^G A i<F , 

2cT^ =Vc (^p A + m A F - (p A b + m A e) A 

2cTj^ = - (i^F) A (m A \/) - {i^ F) A *(m A V) 

2 cr^^ = - Vc (P A + m A F) - K V A (p A b + m A e) - A (p A b + m A e) 

The above are vahd for all simple media in arbitrary gravitational fields. 
For a simple medium, which may be inhomogeneous, anisotropic and intrin- 
sically magneto-electric, at rest in an inertial frame in Minkowski spacetime 
with Minkowski coordinates {t, x} one has (in Euclidean notation) 

V = —dt , e = E ■ dx , h = B ■ dx , h = H ■ dx , d = D ■ dx 
and 

g{e,d)=E-D , g{h,h) = H ■ B 



S = -{E X H) ■ dx 
The coordinate components of the stress-energy-momentum tensor follow as 

Too = |(F ■ D + # ■ 5) 

= -\{EiDj + EjDi) - \{H,B, + B,Hi) + \5i,{E ■ D + H ■ B) (103) 
Tqu = Tko = — (F X H)k 
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11 Conclusions 



Natural assumptions made above for the dependence of the constitutive ten- 
sor Z on the normahsed 4-velocity of a simple medium have led via a non- 
trivial variational argument to a contribution to the stress-energy-momentum 
tensor (involving phenomenological electromagnetic interactions with bulk 
matter) that coincides with that suggested by Abraham under more restricted 
circumstances. Although natural, the assumptions based on physical consid- 
erations are not, however, necessarily the simplest to make. 

If Z is chosen to be independent of the metric and hence V ^ with Zt = 
Zq and Z = so that G = Zq{F) in all gravitational fields, one obtains 
immediately from the above variational calculations (1721) the drive forms 

CTa = 1{F Alai<G -iaG Ai<F) (104) 

and the associated stress-energy-momentum tensor 

T = liaG ® + \iaF ^i^G-^i^iF Ai.G)g (105) 

showing clearly its independence of the 4-velocity of the medium. It is of 
interest to note that such a tensor coincides with that obtained by sym- 
metrising the one proposed by Minkowski. 

In the absence of a generally accepted relativistic covariant description 
of deformable matter interacting with electromagnetic fields, the adoption 
of a particular stress-energy-momentum tensor for the electromagnetic field 
alone in polarisable (and possible magneto-electric) media must remain a 
matter of expediency. However, useful models for the total stress-energy- 
momentum tensor for such systems can benefit from the use of sufficiently 
general phenomenological descriptions of the electromagnetic properties of 
moving media compatible with relativistic covariance. For example a ther- 
modynamically inert (pressureless, cold) fluid can be modelled by adding the 
electromagnetic stress-energy-momentum tensor fl96l) to the matter stress- 
energy-momentum tensor ^ AfV ® V where A/" is a scalar number density 
fleld, rriQ some constant with the dimensions of mass and V the unit time- 
like 4-velocity fleld of the fluid. Supplemented with continuity conditions, the 
vanishing divergence of such a combination yields the dynamics of the sys- 
tem and with prescribed boundary conditions at an interface separating such 
media with different properties one may compute bulk forces and torques. 

A review has also been given of the symmetry constraints expected of the 
total stress-energy-momentum tensor particularly when this is considered to 
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be a source of relativistic gravitation. This led to a definition in terms of a 
variational derivative and a consideration of the response of the electromag- 
netic constitutive properties to gravitational perturbations. It is suggested 
that stress-energy-momentum tensors parameterised by a self-adjoint consti- 
tutive tensor Z offer a viable means to explore the electromagnetic properties 
of a range of inhomogeneous, anisotropic and possibly magneto-electric con- 
tinua, at least in regions where dispersion and losses can be ignored to a first 
approximation. This formulation suggests a method to determine the prop- 
erties of Z by exploring its phenomenological response to electromagnetic 
fields in arbitrarily moving reference frames and variable gravitational fields. 
It opens up the possibility of performing such experiments in new environ- 
ments such as those carried out under terrestrial free-fall or space station 
situations or in astrophysical contexts. 
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12 Appendix 



Using the notation established in the text, this appendix derives the useful 
formula (l63|l relating to T*r\l/ where \l/ G TAPM. Let / denote a multi- 

index constructed from the single indices a, 6, c ... in the range 0, 1, 2, 3 where 
the components of the metric tensor g and ^ in an gf-orthonormal basis {e'^} 
are respectively rjah and \E'/. Thus 

v]/ =\l/^e^ 

and 

ii =i! I + m I {e^y (loe) 

Since is the exterior product of p 1-forms 
Similarly, since the basis is orthonormal 
Thus using (11061) 

<^je^ = - (f M^^ 

Applying * to this gives 

^ = - *(e' A Zc^) (107) 

But 

(^^)- =(^^^eO' 

=^/*e^ + ^/(*e^)' 

Substituting from (I107p yields the relation 

(★^)- = e'' A ic(**) - *(e'' A ic^) + 
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